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We construct the low energy theory of a doped Mott insulator, such as the high-temperature 
superconductors, by explicitly integrating over the degrees of freedom far away from the chemical 
potential. For either hole or electron doping, a charge 2e bosonic field emerges at low energy. The 
charge 2e boson mediates dynamical spectral weight transfer across the Mott gap and creates a 
new charge e excitation by binding a hole. The result is a bifurcation of the electron dispersion 
below the chemical potential as observed recently in angle-resolved photoemission on Pb-doped 
Bi 2 Sr 2 CaCu208 + 5 (Pb2212). 



Two problems beset the construction of a proper low- 
energy theory (explicit integration of the high energy 
scale) of doped Mott insulators. First, the high energy 
degrees of freedom are neither fcrmionic nor bosonic. To 
illustrate, in a Mott insulator, the chemical potential lies 
in a charge gap between two bands that represent elec- 
tron motion on empty (lower Hubbard band, LHB for 
short) and singly occupied sites (upper Hubbard band, 
hereafter UHB). Since the latter involves double occu- 
pancy, the gap between the bands is set by the on-site 
repulsion energy, U. Nonetheless, both double occupancy 
and double holes represent high energy excitations in the 
half-filled insulating state as each is equally far from the 
chemical potential. As neither of these is fcrmionic, stan- 
dard fermionic path integral procedures are of no use. 

Second, unlike the static bands in band insulators, the 
UHB and LHB are not rigid, thereby permitting spectral 
weight transfer. When x holes are placed in a Mott in- 
sulator, at least 2x0| single particle addition states are 
created just above the chemical potential. The devia- 
tion from would be the case in a band insulator, 
is intrinsic to the strong correlations that mediate the 
Mott insulating state in a half-filled band, thereby dis- 
tinguishing Mottness from ordering. The states in ex- 
cess of x arise from two distinct effects. Each hole re- 
duces the number of ways of creating a doubly occu- 
pied site by one, thereby reducing the spectral weight 
at high energy. As the x empty sites can be occupied 
by either spin up or spin down electrons, the 2x sum 
rule is exactpj in the atomic limit. Further, in the pres- 
ence of hybridization (with matrix element t), virtual ex- 
citations between the LHB and UHB increase the loss 
of spectral weight at high energy thereby leading to a 
faster than 2x growth[l], 12, 01 of the low-energy spectral 
weight, a phenomenon confirmed 0, [H, @| widely in the 
high-temperature copper-oxide superconductors. 

Because some of the low-energy degrees of freedom 
of doped Mott insulators derive from the high energy 
scale, low-energy descriptions must cither CI) aban- 
don Fermi statistics or C2) generate new degrees of 
freedom[fj which ultimately leads to electron number 
non-conservation. Current proposals for the low-energy 
physics of doped Mott insulators are based either on per- 



turbation theoryQ followed by projecting out the high- 
energy sector or slaved [8[ operators designed to exclude 
double occupancy. As projection is not integration, nei- 
ther permits an explicit integration of the high energy 
scale, and both miss relevant physical aspects. 

We show that exact integration of the high energy scale 
results in a low-energy theory that possesses a charge 2e 
bosonic mode. Such an excitation might have been antic- 
ipated in light of the mixing between high and low energy 
multiply charged states Our theory is an explicit ex- 
ample of C2 as the conserved charge involves both the 
boson and electron number. Note that the emergence 
of new degrees of freedom in a low energy theory, not 
directly built out of elementary excitations, is not with- 
out precedent. Indeed we believe that there are useful 
lessons to draw for Mott insulators from analogies with 
confining theories or other strongly coupled theories. A 
simple theoretical model which bears some resemblance 
to the theory that we develop below is the non-linear a- 
model (NLcr-M) , in which an initially non-dynamical field 
develops correlations and in fact determines the phase 
structure of the theory. 

While our starting point is the one-band Hubbard 
model, 

#Hubb = -t ]T 9ijc{a c j,<y + U J2 c i,T c <4 c< >i Ci 't W 

our scheme is completely general and is applicable to the 
n-band case as well. Here i,j label lattice sites, <7y is 
equal to one iff i,j are nearest neighbours and Ci a an- 
nihilates an electron with spin a on lattice site i. The 
Hilbert space of this model is a product of Fock spaces, 
(£>i (J-j ® .FjJ. We are concerned in the limit when the 
Hubbard bands are well-separated, U ^> t. Given that 
the chemical potential lies in the gap between such well 
separated bands at half-filling, the high energy degree of 
freedom is ambiguous at half-filling. Both double occu- 
pancy (UHB) and double holes (LHB) are equally costly. 
Doping removes this ambiguity. Hole-doping jumps the 
chemical potential to the top of the LHB thereby defining 
double occupancy to be the high energy scale. For elec- 
tron doping, the chemical potential lies at the bottom of 
the upper Hubbard band and it is the physics associated 
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FIG. 1: Extended Hilbert space (left) which allows an explicit 
integration of the high energy scale. Hopping processes (right) 
included in the Lagrangian. 



with double holes in the lower Hubbard band that should 
be coarse-grained. A generalized particle-hole transfor- 
mation relates the two theories. 

Hole Doping: The basic idea of our construction is to 
rewrite the Hubbard model in such a way as to isolate 
the high energy degrees of freedom so that they can be 
simply integrated out. We do this by first introducing a 
new oscillator that represents the degrees of freedom at 
high energy and including a constraint which ensures that 
the extended theory is equivalent to the Hubbard model. 
If we simply solve this constraint we return to the de- 
scription (1) of the Hubbard model, while if instead, we 
integrate out the high energy degrees of freedom, we will 
obtain the low energy effective theory. To this end, we 
extend the Hilbert space ®i ® -Fj. ® Fd)- We asso- 
ciate D< with the creation of double-occupation. In order 
to limit the Hilbert space to single occupation in the D 
sector, we will take D to be fermionic. Integrating over 
the high-energy scale will be accomplished by integrating 
over D. In particular, we will formulate a Hamiltonian 
for the extended theory in such a way that if we were 
to solve the constraint, precisely the Hubbard Hamilto- 
nian (1) would be recovered. The action of the standard 
electron creation operator, c\ G and the new fermionic op- 
erator, D\ to create the allowed states on a single site 
are shown in Fig. [TJ We now formulate a Lagrangian 
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in the extended Hilbert space in such a way as to in- 
clude the hopping terms (see right panel in Fig. [T]) 
present in the Hubbard model (where we have replaced 
doubly occupied sites by D-occupation). Here, i] is a 
formal complex Grassmann constant which we have in- 



serted in order to keep track of statistics, and cPr] denotes 
Grassmann integration. The parameter V a has values 
V| = 1, Vj. = —1, and simply ensures that D couples 
to the spin singlet. The operator Cij a is of the form 
Cij a = fjrjaija- = fjri{l — ni t v){l—nj t a) with number oper- 
ators nj )CT = cJg.Ci.o-. Note that the dynamical terms that 
appear in the Lagrangian are non-traditional because the 
dynamics with the ci a operators must exclude those sites 
which contain the occupancy cj^c]^. |0). Finally, the con- 
straint H cn „ is taken to be 
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where tpj is a charge 2e bosonic field. The constant s 
will be determined shortly. To see how this constraint 
removes unphysical states that arise from the extended 
Hilbert space, we compute the partition function, 



Z = J[Vc Vc) VD PD f V<p 2V] exp /a Ld \ 



(4) 



in Euclidean signature. The integration over ipi yields a 
series of ^-functions which makes the integral over D triv- 
ial. The resultant Lagrangian given by J d rj W-^Hubb = 
Sicr C J CT c-icr + i?Hubb is identical to that of the Hubbard 
model. This constitutes the ultra-violet (UV) limit 
of our theory. As is evident, in this limit the extended 
Hilbert space contracts, unphysical states such as |1, 0, 1), 
|0, 1,1), 1 1, 1, 1) are set to zero, and we identify |1,1,0) 
with 1 0, 0, 1) . Note there is no contradiction between 
treating D as fermionic and the constraint in Eq. ([3]). 
The constraint never governs the commutation relation 
for D. The value of D is determined by Eq. ([3]) only 
when ip is integrated over. This is followed immediately 
by an integration over D at which point D is eliminated 
from the theory. 

The theory given above permits us to coarse grain the 
system cleanly for U t. The energy scale associated 
with D is the large on-site energy U. Hence, it makes 
sense, instead of solving the constraint, to integrate out 
D. This will result in the low-energy (IR) theory. 
Such an integration may be done exactly as the theory is 
Gaussian in D. This is not possible in previous theories. 
Because of the dynamical term in the action, integration 
over D will yield a theory that is frequency dependent. 
We identify the corresponding low-energy theory by set- 
ting the frequency to zero. Since the theory is Gaussian, 
it suffices to complete the square in the D-field. To ac- 
complish this, we define the matrix 
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and bi = J2j b v — ^2ja Sy Cj.o-Vo-c^ . At zero frequency 
the Hamiltonian [lCj is 

H h R = ~t Yl 9ijOi^cl a c jttT + H int - ^TrlnM 
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which constitutes the true (IR) limit as the high-energy 
scale has been removed. The energy scale s is set by 
noting that the fourth term entering our Hamiltonian 
can mediate spin exchange. As the energy scale for 
this process is t 2 /U , we make the identification s ~ t. 
Hence, appearing at low energy is a charge 2e bosonic 
field which can either annihilate/create doubly occupied 
sites or nearest-neighbour singlets. That the energy cost 
for double occupancy in the IR is t 2 /U and not U under- 
scores the fact that the UHB and LHB are not orthogo- 
nal. If they were, integrating out the high energy scale 
would not result in new charge 2e degrees of freedom at 
low energy. While electron number conservation is bro- 
ken in the IR, a conserved low-energy charge does exist, 
however[llj: Q = J2 ia C L C ^ + 2 X« ViPi- As © im- 
plies, the bosons acquire dynamics only through electron 
motion. Further, they lack a Fock space of their own since 
all operators in the extended space have been integrated 
out. Indeed, on phcnomcnological grounds, weakly inter- 
acting bose-fermi models with a boson Fock space have 
been advanced 12 , [l]| . What the current analysis lays 
plain is that there is a rigorous connection between a 
strongly coupled bose-fermi model and the low-energy 
physics of doped Mott insulators. As we will see, ip^s 
role is to provide internal structure to the electron by 
mediating composite excitations. 

In the limit U — > oo, the theory reduces to the re- 
stricted hopping term and the third term in Eq. 
Performing the ip integration in the partition function, 
we arrive at the constraint 8{c^Cii). This leads to a 
vanishing of double occupancy, the correct result for 
U = oo. Second, for tp — 0, we have the restricted hop- 
ping term and second term in Eq. ©. Approximating 
Mij by its leading term, &y , we reduce the second term to 

EA tfe ; = Ezjw 9ij9ji4,a V ^ c i^' V ^ c ^' which con " 
tains the spin-spin interaction — (Si ■ Sj — TiiUj/A) as well 

as the three-site hopping term. Hence, the p — limit 
contains the t — J model, thereby establishing that the 
physics contained in (pi is non-projective. 
Electron Doping: For electron doping, the chemical 
potential jumps to the bottom of the UHB and hence 
the degrees of freedom that lie far away from the chem- 
ical potential no longer correspond to double occupancy 
but rather double holes. We proceed as before by extend- 
ing the Hilbert space and constructing a new Lagrangian 
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that preserves the distinct hops in the Hubbard model 
where the operator Di is a fermion associated with double 
holes. In this case, the constraint is given by H con — 
sfj^2 i <pi(Di — ?ycj|cjj) + h.c. Two differences to note 
are that 1) because the chemical potential resides in the 
UHB, the electron hopping term now involves sites that 
are at least singly occupied and 2) the order of the Di 
and Ci operators is important. If we integrate over tpi 
and then Di, all the unphysical states are removed and 
we obtain as before precisely LHubb- Hence, both theories 
yield the Hubbard model in their UV limits. They differ, 
however, in the IR as can be seen by performing the 
integration over Di. The corresponding integral is again 
Gaussian and yields 
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as the IR limit of the electron-doped theory. In ad- 
dition to the hopping term, the last term also differs 
from the hole-doped theory as it enters with the oppo- 
site sign. The generalised particle-hole transformation 
(GPHT) that leads to the hole-doped theory is Ci„ — > 
e iQ-r; c t^ (-q _ augmented with ipi — > — p\. As tpi 

is a complex field, the GPHT interchanges the creation 
operators of opposite charge. We again make the iden- 
tification s ~ t because the last term can also mediate 
spin exchange. When the boson vanishes, we do recover 
the exact particle-hole symmetric analogue of the hole- 
doped theory. Because the field tp now couples to dou- 
ble holes, the relevant creation operator has charge — 2e 
and the conserved charge 11 is Q = ^ icr c\ a Ci a — 2<p\<pi. 



This sign change in the conserved charge will manifest 
itself as a sign change in the chemical potential as long 
as (iplcpi) 7^ 0. Likewise, the correct U — » oo limit is 
obtained as before. 

To uncloak how the spectrum changes upon single- 
electron addition or subtraction, we derive an exact ex- 
pression for the electron operator in the new low en- 
ergy theory. To this end, we translate the Lagrangian 
for the hole-doped theory by a source term, that gen- 
erates the canonical electron operator when the con- 
straint is solved. The appropriate source term that 
yields the canonical electron operator in the UV, is 




J2i<r J i," w(l ~ n i,v)4a + VaDjci^T) + h.c. However, 
in the IR in which we only integrate over the heavy degree 
of freedom, Di, the electron creation operator 
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contains the standard term for motion in the LHB, (f — 
n i,a)c\ a with a renormalization from spin fluctuations 
(second term) and a new charge e excitation, c,. ff TW^ 1 ^, 

the IR analogue of the UHB excitation n^g-cj^. Conse- 
quently, we predict that an electron at low energies is in a 
coherent superposition of the standard LHB state (mod- 
ified with spin fluctuations) and a new charge e state 
described by a^M.^ ip'j. It is the presence of these two 
distinct excitations that preserves the 2x sum rule [l| . 

To illustrate this physics, we offer an approxi- 
mate calculation of the electron spectral function. 
For the sake of the following discussion, we con- 
sider ifi to be spatially independent and thereby 
compute the electron Green function by evaluating 
/ d0*d0Fr(/[D<Dci]ci(t)c,-(O)* exp(- / L m dt))/Z. 
Further, we ignore the four-fermion term, b\bi, as this 
term simply renormalizes the standard LHB band as 
Eq. ([9]) indicates. The key results summarised in Fig. 
(J2j) are: I) Below the chemical potential the electron 
spectral function consists of two branches. The inner 
curve corresponds to the standard LHB of the t — J 
model while the outer curve the new excitation arising 
from the non-projective physics in the true low energy 
theory. The bifurcation persists for a wide range of 
doping (1.3 > n > 0.7) and cannot be captured by 
mean-field or saddle-point approximations to the Green 
function. 2) Bifurcation occurs at the second higher- 
energy kink (roughly 0.5i approximately 250meV for the 
cuprates). 3) The difference between the two branches 
is largest at the momentum (0,0) and scales as t. 4) 
The presence of the two branches opens a pseudogap 
in the spectrum as shown in the density of states 
(panel c). The bifurcation, intensity of each branch 
and energy of the kink are in excellent agreement with 
recent experiment 14j. However, the bifurcation 14| has 
been interpreted as evidence for spin-charge separation. 
Although our method is approximate, it is sufficient to 
capture the essence of Eq. ©: two distinct excitations 
(as seen experimentally) constitute the removal of an 
electron at low energies in a doped Mott insulator. 

Whether or not a^MZ 1 ^ constitutes a true bound 
state remains a conjecture at this point, though the pres- 
ence of two states (Fig. in the excitation spectrum 
and Eq. ([9]) are consistent with such physics. An analysis 
based on Bethe-Saltpeter equations is necessary. Conse- 
quently, the physics of doped Mott insulators turns on 
precisely the same kinds of problems that arise in other 
instances of strong coupling such as nuclear structure and 
confinement. 
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FIG. 2: (a) Spectral function along the nodal direction for 
filling n = 0.9 and U = Wt at T = 0. The lower energy 
branch below the chemical potential arises from a binding 
(which opens a pseudogap in the density of states (c)) of a 
hole with the charge 2e boson, (b) Two kinks occur, one at 
0.15t ~ 70meV and a higher-energy kink at .5t ~ 250meV at 
which the bifurcation obtains. 
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